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The Millennium Galaxy Catalogue: Science highlights 
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Abstract. The Millennium Galaxy Catalogue (MGC) provides a structural 
database comprising 10,095 well resolved galaxies drawn from a 37.5 sq. deg re- 
gion of sky, with B MGC < 20.0 mag and 96.1 per cent spectroscopic completeness. 
The data are being used to investigate a number of diverse topics including: the 
nearby galaxy merger rate (via close pairs and asymmetry); dust attenuation; 
bulge and disc luminosity functions; the luminosity size relations; the supermas- 
sive black hole mass function; galaxy bimodality; the space-density of high and 
low surface brightness galaxies; blue spheroids; and the total B-band luminosity 
function. Plans to extend the MGC in area (200 sq. deg), depth (Ky ega = 16.5 
mag), resolution (0.5") and wavelength (u — K) are underway. 



Introduction 



The Millennium Galaxy Catalogue (MGC: iLiske et al.ll2003 ) originated with a 



major award of time with the Wide Field Camera at the Isaac Newton Tele- 
scope. Over a period of three annual cycles (1999-2001) an equatorial strip 75 
deg long by 0.5 deg wide was observed via sequential 12.5 min l?-band expo- 
sures with large overlap regions suitable for bootstrap calibration. All data were 
reduced through the Cambridge Survey Unit Pipeline. Bas ic photomet r y wa s 
performed using SExtractor and calibrated as described in ILiske et al.l ( 20031 ) 



providing photometric and astrometric accuracies of ±0.023 mag and ±0.08" re- 
spectively. The catalogue contains over one million galaxies with 10,095 brighter 
than B UGC = 20 mag over an effective useable area of 30.88 sq. deg, defining 
mgc-bright. All detections were verified by eye and corrected whe re neces- 
sary. The MGC overlaps with th e 2dFGRS NGP (IColless et alJ l20ulh and the 
SDSS Early Data Release region ( Stoughton et al.ll2002l ) providing redshifts for 



almost half of mgc-brigh t, the remainder w ere obtained over the subsequent 
three years (2001-2004, see [5 river et aT1l2005l ) via (predominantly) general user 
2dF time at the AAT, supplemented with additional time at the ANU 2.3m, 
and the NTT, Gemini and TNG telescopes to pursue very high and very low 
surface brightness systems. Most recently we have completed 2D Sersic and 
Sersi c bulge plus exp onential disc profilin g of all 10,095 galaxies (see lAllen et al 



120061 ) with GIM2D (jSimard et aljl2003f h All data products were publicly re 
leased at the IAU General Assembly in Prague and are available via our website: 
http:/ / www.eso.org/~jliske/mgc 
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2. Science highlights 



2.1. The bivariate brightness distribution (M— < fi e ff >-plane) 

In iDriver etHI (12005^ we derive the space-density of galaxies as a function of 
magnitude and mean effective surface brightness (the so called bivariate bright- 
ness distribution). Our 2D variant of step- wise maximum likelihood enables us 
to compensate for the different volumes over which galaxies of similar luminos- 
ity but distinct surface brightnesses can be seen, as well as define the MGC 
selection boundary. The distribution (see Fig. 1) shows the strong luminosity- 
surface brightness relation but more importantly that; (i) galaxies near L* 
(Mb ~ —19.6) are well defined by a Gaussian distribution in surface bright- 
ness which is unimpeded by the selection boundary, and (ii) that faintward of 
Mb = — 16 mag the surface brightness and size selection boundaries start to 
severely hinder our ability to sample the dwarf population, causi ng a loss of 
~_6% of the dwarfs at the high surface brightness end alone (see iLiske et al.l 
2006). The conclusion is that while the bright-end of the galaxy luminosity 
function is well defined, the dwarf regime is poorly constrained and warrants 
further investigation via a deeper survey (both in terms of flux limit, spatial 
resolution and spectroscopic follow up limit). 





M HuD „-5logh (meg) 



Figure 1. The bivariate brightness distribution as a 3D plot (left) and a 
contour plot (right). The solid line denotes the MGC selection boundary. 
Vertical scaling/contours are logarithmic. 



2.2. Galaxy bimodality in colour and structure 

In IDriver etahl (|2006l ) we identified a remarkably strong bimodal distribution 
in the colour and structure plane (see Fig. 2). Subdividing by morphological 
type we found that ellipticals/lenticulars (bulge-dominated systems) lie in the 
red compact peak, late-types (disc only systems) lie in the blue diffuse peak, 
and mid-type spirals (bulge plus disc systems) straddle the two peaks with no 
obvious bimodality. We interpreted this to imply that galaxy bimodality is a 
fundamental reflection of the two-component nature of galaxies (i.e., combina- 
tions of spheroids and discs and not two distinct galaxy types). This in turn 
implies that galaxies may form via a two stage process: spheroids form early 
via a dissipational collapse/rapid merger phase (suggested by their colours and 
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concentrated 3D structure), and galaxian discs grow more slowly via splashback 
(dissipationless collapse), infall and/or accretion (suggested by their colours and 
more fragile 2D structure). In Driver et al. (2007, in prep but see Driver, Liske 
& Graham 2006) we confirm, via our bulge plus disc structural decompositions 
of the mid-type spirals, that the two-component interpretation appears to hold, 
with the bulges and discs of mid-type spirals lying in distinct non-overlapping 
peaks in the colour-structure plane. The current dominant semi-quiescent phase 
is also suggested by our relati yely low merger rate fo und from the MGC data by 
both dynamically close pairs ( De Propris et al. 20051 ) and asymmetry indicators 



(|De Propris et al.ll2007l l. 




Figure 2. The bimodality of galaxies in the colour-structure plane, defined 
by observ ed number, (u — r) C q re an d log(n) (where n is the single profile Sersic 
index, see lGraham fc Driver! 120051 .). The smaller panels show the same data 
but with the z-axis changed to either the volume corrected space-density (left) 
or the stellar mass density (right). 



2.3. The luminosity function of bulges and discs etc. 

In Driver et al.1 ( 2007al ) we derive the luminosity function of discs, red ellipticals, 



classical bulges, blue ellipticals and pseudo-bulges (see Fig. 3). We derive their 
luminosity densities and stellar mass densities and find the following breakdown 
in terms of the stellar mass of 58:13:26:1.5:1.5 (%) respectively. Adopting the 
two-mode model for galaxy formation (rapid bulge assembly followed by semi- 
quiescent disc formation) this implies that the later disc-mode is responsible for 
the bulk of stars visible today (i.e., ~ 60%). 

2.4. Dust attenuation of bulges and discs 

In Driver et al. ( 2007bl ) we derive the bulge and disc luminosity functions in in- 



tervals of inclination and find a strong dependency of the characteristic turnover 
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M B>cc -51ogh (mag) 

Figure 3. The luminosity function of spheroids (red ellipticals+red bulges), 
discs and blue spheroids (blue ellipticals plus pseudo-bulges). 



luminosity with disc inclination. We int erpret this a s due to dust and find our 
data agree well with the predictions of iTuffs et al.1 (|2004l ). Using the Tuffs et 
al. models we constrain the mean face-on opacity of discs and derive a face-on 
attenuation of 0.2 mag for galaxy discs and 0.84 mag for galaxy bulges. To- 
gether the total galaxy attenuation is dramatic (see Fig. 4) implying that only 
50% of the -B-band photons produced by stars emerge from the nearby galaxy 
population with the remainder absorbed by dust grains in the ISM. This has 
dramatic implications for the -B-band luminosity density and derived -B-band 
galaxy luminosity functions, and also affects stellar masses based on optical 
colours. 
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Figure 4. The MGC B UGC band attenuation of discs (left) and bulges (right). 
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2.5. The i?-band luminosity function corrected for dust attenuation 

In Driver et al. (2007, in prep.) we rederive the -B-band galaxy luminosity 
function(s) incorporating our dust attenuation correction (as well as upgrade 
our photometry from Kron to GIM2D magnitudes) and apply individual k(z) 
and e(z) corrections for the bulge and disc components (see Fig. 5). We find 
that the 5-band galaxy luminosity function turnover moves brightwards by 1.0 
mag when dust is correctly taken into account. Returning to the component 
stellar mass contributions we now find that the stellar mass locked up in bulges 
is significantly higher and revise our earlier values to: 58:10:29:1.5:1.5 (see §2.3). 
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Figure 5. The progression of the MGC i?-band luminosity function. 



2.6. Supermassive black holes 

Graham k, Driver] ( 20071 ) provide an updated M^-n relation (akin to the Mbh~ 



CT relation, see Novak, Faber & Dekel 2006) which can be used for predicting 
accurate supermassive black hole masses in other galaxies. Using the measured 
(spheroid) Sersic indices drawn from the 10k galaxies in the MG C, we provide a 



new estimate of the local supermassive black hole mass function (j Graham et al 



2007; see Fig. 6). The observational simplicity of our approach, and the direct 
measurements of the black hole predictor quantity, i.e. the Sersic index, for 
both elliptical galaxies and the bulges of disc galaxies makes it straightforward 
to estimate accurate black hole masses in early- and late-type galaxies alike. 
Moreover, the M^-n relation has a total absolute scatter of only 0.31 dex. 



3. Summary 

The MGC is a relatively modest project but is now starting to produce some 
fairly unique science. In cosmological terms the survey is currently too small to 
place constraints on cosmological parameters but we can investigate the break- 
down of the cosmic baryon budget today. From the MGC we are finding that 
8.3/iq.7% of the baryons produced in the Big Bang are currently in stars, O.OO6/10.7 



6 




lo SlO M BH ( M ©) lo SlO M BH ( M o) lo SlO M BH ( M ©) 



Figure 6. The derived MGC supermassivc black hole mass functions for all 
galaxy types (left), early-types (middle) and late-types (lower). 

per cent are in interstellar dust grains, and 0.008/io.7% are locked up in super- 
massive black holes. Plans are afoot to now expand the MGC by ingesting the 
UKIRT/UKIDSS, VST KIDS and VISTA VIKING data along with extended 
spectroscopic coverage with AAOmega at the AAT and supporting GALEX and 
Spitzer observations. The extended survey (GAMA), will provide a 10 x gain 
in sky coverage, 2x gain in spatial resolution, 2 mag increase in spectroscopic 
follow-up, and 4x gain in wavelength coverage and should provide the ultimate 
post-SDSS structural and spectroscopic galaxy resource for multipurpose galaxy 
studies. One key objective will be to explore CDM on 1 Mpc to 1 kpc scales 
where the interplay between CDM and baryon physics becomes crucial. 
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